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ABSTRACT 

Analytical relations are derived for the amplitude of astrometric, photometric 
and radial velocity perturbations caused by a single rotating spot. The relative 
power of the star spot jitter is estimated and compared with the available data 
for k 1 Ceti and HD 166435, as well as with numerical simulations for k 1 Ceti 
and the Sun. A Sun-like star inclined at i — 90° at 10 pc is predicted to have 
a RMS jitter of 0.087 /xas in its astrometric position along the equator, and 
0.38 m s _1 in radial velocities. If the presence of spots due to stellar activity is 
the ultimate limiting factor for planet detection, the sensitivity of SIM Lite to 
Earth-like planets in habitable zones is about an order of magnitude higher that 
the sensitivity of prospective ultra-precise radial velocity observations of nearby 
stars. 

Subject headings: stars: spots — techniques: interferometric — techniques: ra- 
dial velocities — techniques: photometric — planetary systems — stars: indi- 
vidual (HD 166435, k 1 Ceti) 
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Introduction 



With the anticipated launch of the SIM Lite mission in the near future, we are embarking 
on a long and exciting journey of exoplanet detection by astrometric means. One of the 
main goals of this missi on is the detection of habitable Earth-like planets around nearby 
stars (jUnwin et al.ll2008l ). To date, most of the exoplanet discoveries have been made by the 
Doppler-shift technique, while t he astrometric metho d has been limited to the use of the FGS 
on the Hubble Space Telescope ([Benedict et al.ll2006l ) and to groun d-based CCD observation s 
of low-mass stars with giant, super- Jupiter, planetary companions (IPravdo fc Shaklanll2009l ). 
In achieving the strategic goal of confident detection of rocky, Earth-sized planets in the 
habitable zone, the prospective astrometric and spectroscopic ultra-precise measurements 
will encounter a number of limitations of technical and astrophysical nature. 

For the Doppler-shift technique, many of these limitations will be dealt wit h by further 
impr ovement in the instrumentation or refinement of the observational procedure (IMayor fc Udry 
20081 ). However, the presence of astrophysical noise due to stellar magnetic activity emerges 
as the ultimate bound on the sensitivity of planet detection techniques, and the only remedy 
suggested thus far, is selection of particularly inactive, slowly rotating stars. Indeed, a very 
small fraction of stars in the high-precision HARPS program of exoplanet search exhibit ra- 
dial velocity (RV) scatter of less than 0.5 m s -1 . Although this type of variability is probably 
driven by the rotation of bright and dark structures on the surface (star spots and plage 
areas), the frequency power spectrum of such perturbations can b e fairly flat, exten d ing to 
frequencies much higher or lower than the rotation, as shown by ICatanzarite et al.l (120081 ) 



for the Sun. Arguments have been presented (e.g., Eriksson fc Lindegrenl 120071 ) that star 
spots also result in very large astrometric noise of ~ 10 /zAU, which should thwart discovery 
of habitable Earth analogs. The aim of this paper is to quantify the effects of rotating spots 
in astrometric photometric and RV measurements more accurately, taking into account the 
limb darkening, geometric projection and differential rotation, and to assess the expected 
vulnerability of the RV and astrometric methods to such perturbations. We do this by direct 
analysis as well as by numerical simulation, and support our findings by the data for the 
Sun and two rapidly rotating stars. 



2. Perturbations from a single spot 

We consider a single circular spot whose instantaneous position on the surface in the 
stellar reference frame is given by longitude I and latitude b, which are the angles from 
the direction to the observer and from the equator, respectively. The projected area of the 
spot is 7rrcosC, where C is the central angle between the direction to the observer and 
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the center of the spot, and r is the radius of the spot in radians, r << 1. The position 
vector of the spot in the local sky triad {9ft, E, N} is s = s 2 , s 3 ] T = [— cos/ cos b smi — 
sin b cosi, — sin/ cos b, — cos / cos b cosz + sin b sini] T . N is north, E is east, and 9ft is the 
line-of-sight directions in this right-handed triad. The velocity vector of the spot is V = 
[t>i, t>2, t>3] T = [sin / sin i, — cos/, sin / cosi] T Vb, where the differential rotation velocity 

2irR cos b 

Prot{b) 

where V eq is the equatorial rotation velocity. This equat ion assumes the differential relation 
for th e Sun derived from sunspot latitudes and periods (INewton fc Nunnlll95ll ; iKitchatinov 



V b = - w K q cos b (1 - 0.19 sin 2 b), (1) 



20051 ). We also assume that the spot's contrast with respect to the local surface brightness is 



fixed at / s . The asymmetry in the distribution of surface brightness due to the spot results 
in certain perturbations in the integrated flux, photocenter and radial velocity of the stellar 
disk. Assuming the limb darkening for the Sun at A = 550 nm 

1(C) 

JL-Z = o.30 + 0.93 cos C-0.23 cos 2 C, (2) 

the integrated flux from the stellar disk is 0.9057r/(0), the following relations are obtained 
for the amplitudes of perturbation 
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where R is the apparent radius of the star. These expressions describe the modulation of 
the flux, photocenter, and radial velocity of the star due to the motion of a spot. For the 
Sun, = 2 km s _1 and R@ = 4650 fiAU. We estimate a relative flux variability of the Sun 
of RMS(A/ 1 // 1 ) = 3.24 • 10 -4 afte r subtracting a 10-yr p eriod solar cycle light curve from 
the solar irradiance PMOD data (IFrohlich fc Leanl Il998l ). Therefore, the sunspot-related 
jitter is not greater than AmR (1.5 /xAU for the Sun) in position and than AmV eq (0.65 m 
s _1 for the Sun) in radial velocities, where Am is the characteristic magnitude jitter. The 
astrometric perturbation decreases with distance to a typical value of 0.15 /xas for the Sun 
at D — 10 pc. 

Eqs. |3] can be used for numerical simulation of star spot perturbations in a computa- 
tionally efficient way. They can also be integrated in quadratures to estimate the power 
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(root-mean-square, RMS) of the jitter. This results in fairly tedious series in powers of sines 
and cosines of b and i, which we do not give here for brevity. Some of the results are rep- 
resented in graphical form in Fig. [TJ It should be noted that the ratios of RMS values in 
flux and each of the remaining parameters can only be directly computed for a single spot 
on the surface. Using this figure, and the RMS jitter for the Sun at i = 90° and 6 = 0° are 
RMS(Ax) = 0.87 /xAU and RMS(AVr) = 0.38 m s" 1 . 

The geometric projection factors S2 and v\ cos b differ only by a constant factor — sin i. 
Therefore, the ratio of the perturbation amplitudes, as well as perturbation RMS in Ax and 
in AVr is constant for a single spot: 

RMS(AVr) Kqsini niQ . a ., m 
RMS(Axy = ^^ (1 - - 19sm ^ (4) 

Using the values for the Sun, the approximate scaling relation is 
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in units m s _1 /ias -1 , where P ro t,© = 24.47 d is the sidereal rotation period of the Sun at the 
equator. Using the Carrington period of 25.38 d instead will to some extent account for the 
distribution of sunspots in latitude, and allows one to ignore the differential rotation factor 
to first-order approximation. 

The data in Fig. CD and Eq. [5] can be used to estimate the relative magnitude of starspot 
jitter in astrometry and RV measurements. Eqs. [3] provide an efficient and direct way of 
simulating these perturbations for any configuration of spots. 



3. Comparison with observations 
3.1. HD 166435 

The star HD 166435 is a solar-type dwarf at 25 pc without conspicuous signs of chro- 
mospheric or coronal activity, which nonetheless ex hibits large and cor related variations in 



brightness, radial velocity, and Call H and K lines (IQueloz et al.ll200ll ). Strong evidence is 
presented in that paper that these periodic variations are caused by a photospheric spot or a 
group of spots, including analysis of spectral line bisectors. Some confusion with a possible 
short-period giant planet resulted from the conspicuously large amplitude of the RV varia- 
tion (~ 200 m s~ x peak-to-peak) and the stable phase on a time scale of 30 days. Queloz 
et al. obtained a projected rotation velocity of v sin i = 7.6 ± 0.5 km s _1 and a period of 



50 100 150 50 100 150 

inclination (i), degr inclination (i), degr 



Fig. 1. — The magnitudes of relative RMS perturbations from a single star spot: a) ra- 
tio RMS(Ax)/ RMS(AF/F) in units of apparent stellar radius R; b) ratio RMS(AVr)/ 
RMS(AF/F) in units of V eq (l - 0.19 sin 2 b). In both cases the values at % = 90°, 6 = 0° are 
0.448. 
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3.7987 d. They estimated a i = 30° from these data. The shape of the variability curves is 
consistent with a single dominating spot rotating with this period. 

Ignoring the differential rotation term (which is probably small for fast-rotating stars) 
and setting AVr to 200 m s _1 from Queloz et al.'s Fig. 8, or 166 m s -1 from their text, 
and Am to 0.1 mag or 0.07 mag, we derive from Eqs. [3] cos b = 0.286, b = 73°, or cos b = 
0.339, b = 70°, respectively. This latitude is ambiguous with respect to sign, the possible 
combinations being % = 30°, b = 70° (counterclockwise rotation), or i = 150°, b = —70° 
(clockwise rotation). In any case, the center of the spot is close to the pole, and because 
of the small inclination, circles c lose to the m i ddle of the visible stellar disk. This is fully 
consistent with the conclusion by I Queloz et al.l (120011 ). which they draw from the smoothness 
of the variability curves. Finally, using the above estimates for i and b, we compute the light 
curve, which is indeed a smooth sinusoid-like function of time, with a peak-to-peak amplitude 
of 0.51(1 — f s )r 2 . For a large spot area, the factor r 2 is interpreted as the fraction of the 
observed hemisphere covered by the spot. The characteristic contrast ratio of sunspots is 
0.2 in the optical passban d, which correspon ds to an effective temperature of ~ 4400 K for 
the spotted photosphere (lLanza et al.l 120081 1 . This yields a striking number for the area of 



the spot, r 2 = 0.23 , or an angular rad ius p = 40°. This feature appears to be a dark sea 
engulfing the pole. iQueloz et al.l (120011 ) find that the phase of the RV variation is confined 
to a ±0.1 interval (their Fig. 4), which indicates that the feature is stable on a time scale of 
2 years, but with considerable internal variations of brightness. 



3.2. k 1 Ceti 



The available data for this star relevant to this study include the precision photometric 
series from MOST and 44 individual RV measurements spread over some 20 years, provided 
by one of us (DF). The MOST data sets have been carefu lly analyzed in other papers 



(IRucinski et al.ll2004l ; I Walker et al. 1 120071 : iBiazzo et al.ll2007l ). with a firm conclusion that 



its light curve can be well modeled with a set of two or three dark spots. The periods of 
rotation range between 8.3 and 9.3 days. The observed peak-to-peak amplitudes in flux are 
roughly 0.05 in 2003, but only ~ 0.0 2 in 2004 and 2005. The projected speed of rotation 
measured by (jValenti fc Fischerl 120051 ). v sin % = 5.2 ± 0.5 km s _1 implies, from Eq. [3j a RV 
perturbation < 260 m s _1 in 2003 and < 120 m s _1 in 2004 and 2005. Due to the sparsity 
of the RV data, we can not match them directly with the intervals of MOST observations, 
but we can estimate the amplitude of RV variation over two decades. With the smallest 
observed value of —68.56 ± 5.01 m s -1 , and the largest +42.50 ± 9.88 m s _1 , the amplitude 
is close to the single-spot estimate for the more quiescent periods, but is much smaller than 
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the prediction for 2003. Independent observations by I Walker et al. I (119951 ) during 1982- 
1992 also indicated a peak-to-peak amplitude of ~ 100 m s _1 . One possible explanation is 
that the spots usually reside at high latitudes, which reduces the relative variation in RV 
because of the cos b factor. The smaller RV variability may also be related to the fact, that 
two or three spot groups are present on the surface at a time, rather than one single spot. 
Eqs. |3] can not be simply scaled for the case of multiple spots, because the AVr and Ax 
function of time are symmetric around the central meridian (C = 0). As a result, two spots 
well separated in longitude can counterbalance eac h other, reducing or n early canceling the 
net perturbation. Indeed, the detailed modeling by I Walker et al. I (120071 ) indicate that large 
spots frequently occur in the near-polar regions of n 1 Ceti, and that two or three coexisting 
spots are spread in longitude, rather than grouped in a confined active area. 
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Fig. 2. — Numerical simulations of variation in relative flux, equatorial shift of the photo- 
center, and radial velocity of n 1 Ceti caused by two rotating spots, corresponding to the 
first segment of obse rvations with MOST in 2003. Our prediction is drawn with a solid line, 
Walker et al. I (120071 ) results with dashed line, and MOST data with asterisks. 
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Learning more about the properties of spots on this star requires more accurate com- 
putation. One of us (JL) created a code to model differentially rotating star spots and the 
resulting perturbations of the observable parameters by pixelization of the stellar surface and 
integration over the visible hemisphere. The free model paramet ers were optim i zed o n the 
photometric series of k 1 Ceti, essentially repeating the study by I Walker et al. I (120071 ) , but 
also extending it to astrometric and radial velocity predictions. A detailed description of this 
model will be published elsewhere, here we only discuss some of the results relevant for this 
paper. Fig. [2] shows the simulated perturbations from our model, along with the expected 
perturbations from the model by Walker et al., and the actual light curve from MOST for 
the segment of 2003. With only two spots in both cases, the goodness of photometric fit is 
similar with the two models, but our simulation predicts a smaller amplitude of RV variation 
(~ 110 m s" 1 ), probably because of a more symmetric configuration of spots. This prediction 
is in fact consistent with the available RV data. The ratio of simulated variations in AVr 
and Ax is consistent with Eq. [5j 



4. The Sun 

We estimated in §[2] that the expected RMS jitter in RV for the Sun is 0.38 m s _1 . The 
half-amplitude of the reflex motion caused by the Earth orbiting the Sun is slightly less than 
0.1 m s _1 . Does this result imply that detection of Earth-like planets in the habitable zone 
of Sun-like stars is impossible? This brings up the more subtle issue of the spectral power 
distribution of starspot jitter. The characteristic period of star spot rotation is 1 month 
for Solar-type stars, but the orbital periods of habitable planets are of order 1 year. It is 
therefore not obvious without more detailed analysis that the signal-to-noise ratio will be 
too small for a confident detection in the frequency domain of interest. 

One of us (JC) perfo rmed extensive Monte-C arlo simulations for the Sun, which are 



described in more detail in (ICatanzarite et al.ll2008l ). Sunspot groups are generated randomly 



through a Poisson process, with probability distributions consistent with the current data. 
The main purpose of this simulation is to faithfully reproduce the power spectrum of the 
solar irradiance data on a time scale of 30 years. The average number of sunspot groups and 
the dispersion of lifetimes was adjusted in such a way that the predicted and the observed 
spectral power of AF/F match closely in the frequency range 10~ 7 to 4 x 10~ 8 Hz. This gives 
some assurance that the model predictions can be accurately made for the power spectra 
of the astrometric and RV jitters. Typical RMS values of jitter predicted by the numerical 
simulation at i — 90° are: 0.39 m s" 1 for RV, 0.97 fiAXJ in x, and 0.52 /xAU in y, quite close 
to the analytical estimates in § |2j 
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Fig. 3. — Square root of the spectral power of simulated astrometric (left plot) and radial 
velocity (right plot) perturbations of the Sun seen at inclination 90°. The X axis is aligned 
with the equator, and the Y with the rotation axis. 

Fig. [3] shows the square root of power in 100 astrometric and radial velocity observations 
over 5 years simulated for the Sun at i — 90°. The equator is expected to be coplanar with 
the orbital plane, hence, the exoplanet signal will be present only in the x-measurements. 
The power of spot-related jitter spreads far and wide from the rotation period of 25 days. 
It peaks between periods of 1 and 2 years, reaching almost 0.18 /iAU in astrometry. In the 
worst case, exoplanets with signatures of 0.63 /iAU or greater can be confidently discovered 
with SNR> 3.5 by the astrometric method. The spectrum of the simulated RV variations is 
practically identical to the spectrum of x-jitter, as predicted in § [2j The peak value is 0.07 
m s -1 . The corresponding semiamplitude of exoplanet signature detectable at SNR= 3.5 is 
0.25 m b _1 . 



5. Conclusions 

Our results for the Su n are in good agreement with the approximate relations by 



Eriksson fc Lindegrenl (120071 ). who estimated a positional standard deviation of 0.7 /iAU. 
At the same time, their conclusion that "for most spectral types the astrometric jitter is ex- 
pected to be of the order of 10 micro- AU or greater" is misleading, because it is largely based 
on overestimated values of photometric variability from ground-based observations, and it 
does not differentiate the luminosity classes of giants and dwarfs. It can not be concluded 
that the Sun is exceptionally inactive compared to its peers just because the ultra-precise 
solar irradiance data, such as PMOD or SOHO, reveal a much smaller scatter than the infe- 
rior photometric data for other stars. We investigated the indices of chromospheric activity 
(\ogR' UK ) and available rotation periods for some 80 SIM targets, and found that half of 
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them should rotate with the same rate as the Sun, or slower. 



Hall et al.1 (120001 ) presented a detailed study of variability of solar-type stars and its 



relation to the index of chromospheric activity logi?H K , based on 14 years of photometric 
and spectroscopic observations. They found that the Sun at log-R^K = — 4.96 is not more 
variable than its F-G peers at the low end of the activity distributi on. Given that most 
Solar- type stars exhibit similar low levels of chromospheric activity (IGray et al.1 120031 ). we 
expect that finding stars with levels of jitter similar to, or lower than the Sun, should not 
be a problem. Low-jitter, stable stars are common an d plentiful, which a ugurs well for the 
prospects of finding small, rocky planets with Kepler ( IBatalha et al.ll2002l ). 



The SIM Lite Astrometric Observatory (formerly known as the Space Interferome- 
try Mission) will achieve a si ngle- measurement a ccuracy of 1 na . s or b etter in the differ- 
ential regime of observation (jUnwin et al.1 120081 ; iDavidson et al.l |2009i ). Several previous 
studies have addressed SIM' s exoplanet detection and orbital characterization capabilities 
(ICatanzarite et al.ll2006l . and references therein). The "Tier 1" program includes ~ 60 nearby 
stars for which the astrometric signature of a terrestrial habitable planet is lar ge enough to 
be co nfidently measured by SIM. The recently completed double-blind test (ITraub et al. 
20091 ) demonstrated that Earth-like planets around nearby stars can be discovered and mea- 
sured even in complex planetary systems. In this paper, we are concerned with the more 
general question of the ultimate limit to planet detectability set by the activity-related jitter. 
Assuming that the instrumentation progresses to levels where observational noise becomes 
insignificant, which technique holds the best prospects for detection of habitable Earth-like 
planets? Table [Q summarizes the expected RMS jitter and the exoplanet signal for three typ- 
ical nearby stars. In all cases, the solar spot filling factor (r 2 ) is assumed. The signal-to-noise 
ratios (SNR) per observation include only the star spot jitter. Note that the astrometric 
SNR in this case is independent of the distance to the host star, because both the signal 
and the star spot jitter are inversely proportional to distance. As a measure of the relative 
sensitivity of the two methods, the ratio of the SNR values for a given star is independent 



Table 1. Observable signals and star spot jitters for an Earth-like planet orbiting a typical 

dwarf star at 10 pc. 



Star type 


Sun 


F5V 


K5V 


Rotation period, d .... 


25.4 


18 


30 


Astrometric signal, fias 


0.30 


0.23 


0.45 


RV signal, m s~ 1 


0.089 


0.078 


0.109 


Astrometric jitter, fias 


0.087 


0.113 


0.063 


RV jitter, m s - 1 


0.38 


0.69 


0.23 


Astrometric SNR 


3.4 


2.0 


7.1 


RV SNR 


0.23 


0.11 


0.47 
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of the planet mass or the filling factor to first-order approximation. The physical radius of 
the star and the period of rotation are the two parameters with the largest impact on the 
relative sensitivity, but their combined effect is rather modest for normal stars, as is seen 
in the Table. Therefore, in the ultimate limit of exoplanet detection defined by intrinsic 
astrophysical perturbations, the astrometric method is at least an order of magnitude more 
sensitive than the Doppler technique for most nearby solar-type stars. 

The authors thank G. Walker for his detailed and helpful review. The research described 
in this paper was in part carried out at the Jet Propulsion Laboratory, California Institute 
of Technology, under a contract with the National Aeronautics and Space Administration. 
This research has made use of the NStED database, maintained at NExScI, Pasadena, USA. 
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